Vortex Dynamics in Percolative Superconductors 
Containing Fractal Clusters of a Normal Phase 
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The effect of fractal clusters on magnetic and transport properties of percolative superconductors 
is studied. The superconductor contains percolative superconducting cluster carrying a transport 
current and clusters of a normal phase. It is found that normal phase clusters have essential fractal 
features. The fractal dimension of the boundary of normal phase clusters is estimated. The current- 
voltage (V-I) characteristics of superconductors containing fractal clusters are obtained. It is found 
that the fractality of the cluster boundary intensifies pinning. This feature permits to enhance the 
current-carrying capability of the superconductors. 
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I. INTRODUCTION 

An essential feature of the clusters of defects in su- 
perconductor lies in their capability to trap a magnetic 
flux.— i2iMi& The cluster structure affects the vortex dy- 
namics in superconductors, especially when clusters have 
fractal boundary^^ i 10 ' 1 1 A prototype of superconduc- 
tor containing inclusion of a normal phase is a supercon- 
ducting wire. 

The first generation high-temperature superconductor 
(HTS) wires are fabricated following the powder- in-tube 
technique (PIT). The resulting product consists of one 
or more superconducting cores armored by the normal 
metal. The sheath endows the wire with the necessary 
mechanical and electrical properties. The best results 
are obtained now for the silver-sheathed bismuth-based 
compositesj 12 ! 13 In view of the PIT peculiarity the first 
generation HTS wire has a highly inhomogeneous struc- 
ture. Superconducting core represents a dense conglom- 
eration of BSCCO micro-crystallites containing normal 
phase inclusions inside. 

The second generation HTS wires have multi-layered 
structure consisting of the metal substrate (nickel- 
tungsten alloy), the buffer oxide sub-layer, HTS layer 
(YBCO), and the protective cladding made from the no- 
ble metal (silver). In the superconducting layer there are 
clusters of columnar defects that can be created during 
the film growth process as well as by heavy ion bom- 
bardment. Such defects are similar in topology to the 
vortices, therefore they suppress effectively the flux creep 
that allow to get the critical current up to the depairing 
value jii 2 -! 3 - 



II. VORTEX DYNAMICS IN PERCOLATIVE 
SUPERCONDUCTORS 

A superconductor containing isolated clusters of a nor- 
mal phase allows for effective pinning, because the vor- 
tices cannot leave them without crossing the supercon- 
ducting space. The clusters consist of sets of normal 



phase inclusions, united by the common trapped flux 
and surrounded by the superconducting phase. 8,9 When 
the current is increased the vortices start to break away 
from the clusters of pinning force weaker than the Lorentz 
force. Then, the vortices will pass through the weak links, 
connecting the normal phase clusters. In this case depin- 
ning has percolative character ) 14 ' 15 because vortices move 
through randomly generated channels. Weak links form 
readily in HTS's due to the intrinsically short coherence 
length^ 

Let us take the area of cluster cross-section by the 
plane carrying the transport current as a measure of its 
size. Magnetic flux trapped into a cluster is propor- 
tional to its area. Hence the decrease in the trapped 
flux A$ can be expressed with the probability W (A) = 
Pr{VA,- < A} of finding the clusters of area Aj smaller 
than a given magnitude A: 



— = 1-W(A) 



(1) 



According to weak link configuration each nor- 
mal phase cluster has its own value of dcpinning current, 
which contributes to the overall distribution. Thus the 
decrease in the trapped flux is proportional to the num- 
ber of all the normal phase clusters of critical currents 
less than a preset value / and can be expressed with the 
probability F (I) = Pr {Wlj < 1} : 



F(I) 



(2) 



The critical current and cluster area distributions are 
interdependent, because larger cluster has more weak 
links over its boundary and, consequently, the smaller 
dcpinning current. 



III. WEAK LINK DISTRIBUTION OVER THE 
CLUSTER BOUNDARY 

Let us analyze how the vortex exits from a normal 
phase cluster. As the transport current increases, the 
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Lorentz force, pushing the vortex out, will increase as 
well. In order to leave the cluster the vortex has to en- 
ter into one of the weak links, which are randomly ar- 
ranged along the cluster perimeter. Suppose that after 
the vortex reaches the entry point it passes all the way be- 
tween two adjacent normal phase clusters without being 
trapped inside the weak link itself. The vortex exit can be 
considered as the result of random walks under the action 
of the Lorentz force. The following outcomes of the ran- 
dom walks may happen: (a) the vortex enters the weak 
link and leaves the normal phase cluster, (b) the vortex 
does not enter the weak link and continues its random 
walks, and (c) the vortex does not enter the weak link 
at all and remains to be locked inside. The mean num- 
ber of the entry points on the cluster perimeter gives the 
probability measure of the random walk outcomes, which 
are favorable for the vortex to go out. So an exit of the 
vortex from a normal cluster can be treated as the two- 
dimensional generalization of the problem of a random 
walk particle reaching an absorbing border. Unlike the 
classic problem of the distribution of the exit points ^ 
here the boundary of the area is not absorbing all over, 
but there are only discrete absorption points where the 
vortices can enter the weak links. Moreover, the situa- 
tion is complicated by the fact that a random walker is 
permanently subjected to the Lorentz force. 

The distribution of entry points over cluster perimeter 
varies from one cluster to another, so that each cluster 
has the entry point distribution function ip (I) of its own, 
where I is the co-ordinate along perimeter. The probabil- 
ity distribution of functions i/j (I) over all the clusters can 
be characterized by the probability Pr{ip (I)} of finding 
a given function ip (I). 

The most probable function of entry point distribution 
is the mean over all functions 



If all the entry points are uniformly distributed over 
the cluster perimeter, the functional takes the form 



^(0 = J D^(l)^(l)Pr{^(l)} 
(fi) 



(3) 



The path integral Fourier transform on the probability 
functional Pr{ip(l)} represents the characteristic func- 
tional 

/ Vip (I) exp (ifdlk (I) i/> (I)) Pr {4, (I)} 

(O) 

(4) 

where k = k (I) is the element of a reciprocal function 
space. 

If all the clusters are of an equal entry point distri- 
bution, which coincides with the most probable one ([3]), 
the probability Pr{?/>(7)} is zero for all ip (I) that differ 

from ^>(Z), whereas Pr|-0(7)j = 1 . At that rate the 

functional (|4]) becomes 



H[k (1)} =cxp fi^- jdl 



k(l) 



(6) 



where the constant (3 is being chosen to normalize the 
distribution function i\) (I) to unity, so that j3N = 1. The 
functional ((6|) has the form of equation ((5J for the uniform 
distribution of entry points: i/j (I) = This means 

that all the clusters have the same uniform distribution 
of the entry points, for which the probability of finding 
a weak link at any point of the perimeter is independent 
of its position. 

If the concentration of entry points per unit perime- 
ter length n = N/P is constant for all clusters, and all 
the clusters are statistically self-similar, the mean num- 
ber of entry points N along perimeter is proportional to 
its length TV = §n(T)dl = nP . The more entry points 
into weak links are accessible for random walk vortices 
driven by the Lorentz force, the more is the probabil- 
ity that the vortex will leave the cluster, and therefore, 
the smaller is the Lorentz force required to push the vor- 
tex out. Hence, we may write the following relationship 
between the critical current of the cluster, at which the 
magnetic flux ceases to be trapped inside, and its geo- 
metric size I oc 1/N oc 1 jP . 

Thus, to deal with the distribution function ([1]), the 
relation between perimeter and area of clusters should 
be studied. It might be natural to suppose that the 
perimeter-area relation obeys the well known geometric 
formula P oc y/A. However, it would be a very rough 
approximation, because this relationship holds for Eu- 
clidean geometric objects only. As it was first found in, 8 
the normal phase clusters can have fractal boundaries, i. 
e. the perimeter of their cross-section and the enclosed 
area obey the scaling law 



P~b cx A? 



(7) 



H [k (I)} = exp [i 4> dlk {I) V (0 



(5) 



where D is the fractal dimension of the cluster 
boundary^ The fractal nature of such clusters affects 
the vortex transport and depinning in superconductors.- 



IV. FRACTAL GEOMETRY OF NORMAL 
PHASE CLUSTERS 

In order to clear up how the developed approach can 
be used in practice, the electron photomicrographs of 
YBCO films have been studied. The films were pre- 
pared by magnetron sputtering on sapphire substrates 
with a cerium oxide buffer sublayer. The normal phase 
clusters were formed by columnar inclusions of nonstoi- 
chiometric composition. These inclusions were created at 
the sites of defects on the boundary with the substrate 
in such a way that they were oriented normally to the 
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FIG. 1: Perimeter-area relationship for the normal phase clus- 
ters with fractal boundary in YBCO film. The solid line in- 
dicates the least square regression line. 



plete gamma function, A$ and g are the parameters of 
gamma distribution that control the mean area of the 
cluster A = [g + 1) Aq and its variance a 2 A = (g + 1) Aq. 
Exponential distribution is the simplest case of gamma 
distribution of g = 0. 

In accordance with starting formulas (fTJ) , |(3J) as well 
as with scaling law ([7]) we can get the relation between 
the critical current of the cluster and its geometric size: 
/ = aA~ D / 2 , where a is the cluster form factor. The 
cluster area distribution ([5]) gives rise to the critical cur- 
rent distribution 

F( l ) = (T(g + l))- 1 r[g + l 7 G l - 2 / D ) (9) 

where G = {6° / (6 9+1 - (D/2) exp (9) T (g + 1, 8))) 2 ^ D , 
9 = g+l + D/2, r (V, z) is the complementary incomplete 
gamma function, i = I / I c is the dimensionless electric 
current, I c = a (AqG) D ^ 2 is the critical current, which 
gives the point of intersection of the current axis and the 
tangent line drawn through the inflection point of the 
dependence of differential resistance on the current. The 
found distribution allows to derive the probability 
density f(i) = dF/di for the critical currents 



film surface. The profiles of cluster sections by the film 
plane were clearly visible on the photomicrographs, and 
their perimeters and enclosed areas have been measured. 
So the cross-sections of the extended columnar objects, 
which the normal phase clusters are, have been investi- 
gated. The normal phase has occupied 20% of the total 
surface so the percolative superconducting cluster was 
dense enough. The perimeters and areas of clusters have 
been measured by covering their digitized pictures with a 
grid of spacing 60x60nm 2 . The sampling has contained 
528 normal phase clusters located on the scanned area 
of 200^m 2 . A high skewness (1.765) as well as the sta- 
tistically insignificant (5%) difference between the mean 
cluster area and the standard deviation has confirmed 
that the distribution of the cluster areas is exponential. 

The obtained data are given in Fig. Q] All the points 
fall on a straight line on double logarithmic scale with 
correlation coefficient of 0.929. This graph shows that the 
perimeter-area scaling relation, which is inherent to frac- 
tals, is valid in the range of almost three orders of mag- 
nitude in cluster area. The slope of the perimeter-area 
regression line gives the fractal dimension of the cluster 
boundary D = 1.44 ±0.02. This fact that the scaling law 
([7]) with fractional exponent D is fulfilled for the normal 
phase clusters, gives an evidence for their fractal nature. 

In the general way the cluster area statistics may be 
described by gamma distribution 



W(A) = (T(g + l)y 1 1 (g + l 7 ^J (8) 
where T (is) is Euler gamma function, 7 (y) is the incom- 



2G g -* 



DT(g+l) 



r (2/U)( 9 +l)-l exp ^_ Gr 2/^ (1Q) 



V. RESISTIVE STATE OF 
SUPERCONDUCTORS WITH FRACTAL 
CLUSTERS OF A NORMAL PHASE 

Each normal phase cluster contributes to the total crit- 
ical current distribution, so the voltage across a sample 
can be represented as the response to the sum of effects 
from each cluster 



V 



R 



f 



{i - i') f {%') di' 



(11) 



where Rf is the flux flow resistance. Thus, using the 
critical current distribution (fit)]) , we get the V-I charac- 
teristics: 



— = , 1 J ir(g+l,Gi- 2 ' D ^ 
R f T(g + 1)[ \ y J 



G D l 2 T[q + l-^-,Gi- 2 l 



(12) 



The V-I curves for the different fractal dimensions are 
presented in Fig. [21 The lines drawn for Euclidean clus- 
ters (D — 1) and for the clusters of the most fractality 
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FIG. 2: Current-voltage characteristics of superconductor 
containing fractal clusters of a normal phase. The inset shows 
the dependence of differential resistance on a transport cur- 
rent, which is typical for the vortex glass. 



(D = 2) bound the region the V-I characteristics can 
fall within. The curve between them gives the case of 
the previously found fractal dimension D = 1.44. The 
inset in Fig. [2] shows the dependence of differential resis- 
tance Td = dV/di on the transport current. Differential 
resistance is proportional to density of vortices n broken 
away from the pinning centers: = i?/n$o/ B, where 
B is the magnetic field, $0= hc/(2e) is the magnetic 
flux quantum, h is Planck constant, c is the velocity of 
light, and e is the electron charge. It is just a motion of 
free vortices induces electrical field. The dependencies of 
resistance on the current shown in this graph are typi- 
cal for the vortex glass: the curves have a convex form 
on double logarithmic scale as well as resistance tends 
to zero with decreasing current as a result of flux creep 
suppression i 19 i 20 

All the V-I characteristics in Fig. [5] are virtually start- 
ing with the point of i = 1. It is seen that the clus- 
ter fractality reduces electric field arising from the vor- 
tex motion. The reason of this phenomenon lies in the 
peculiarity of the critical current distribution (fT0| . As 
the fractal dimension increases, this distribution broad- 
ens out, moving towards greater magnitudes of current. 
It means that more and more of the small clusters, which 
can trap the vortices best, are being involved in the pro- 
cess. The smaller part of the vortices can move, the 
weaker the induced electric field. Resistive transition 
widens, shifting towards higher currents. The effect of 
fractal dimension can be characterized by the pinning 
gain factor fc* = 20 log (A$ (D = 1) /A$ (£>)), which is 
equal to relative decrease in the number of vortices bro- 
ken away from the clusters of fractal dimension D com- 
pared to the Euclidean ones, as well as by the voltage 



FIG. 3: Pinning gain for the different fractal dimensions of 
the cluster boundary. In the inset the voltage attenuation 
factor is shown. 



attenuation factor k v = 20 log (V (D) /V (D = 1)), char- 
acterizing reduction of electric field. The dependencies of 
these factors on the transport current for different frac- 
tal dimensions are shown in Fig. [3] It is necessary to 
keep in mind that these curves, as well as the graphs 
inFig. [2j are plotted versus dimensionless current nor- 
malized relatively to the critical current J c , the value of 
which depends on the fractal dimension. The pinning 
gain characterizes the properties of a superconductor in 
the range of the transport currents i > 1. At smaller cur- 
rent the breaking of the vortices away has not started yet 
for lack of pinning centers of such small critical currents. 
The pinning enhancement due to the cluster fractality in 
the neighborhood of resistive transition can only be re- 
alized in the case of efficient heat removal that prevents 
the development of thermo-magnetic instability. As for 
any hard superconductor the energy dissipation in the 
resistive state does not mean the destruction of phase 
coherence yet. Some dissipation always accompanies any 
motion of vortices that can happen in even at low trans- 
port current. Superconducting state collapses only when 
a growth of dissipation becomes avalanche- like as a result 
of thermo-magnetic instability. 



VI. CONCLUSION 

In the present work the fractal nature of the normal 
phase clusters is revealed. It is found that the fractality of 
cluster boundary strengthens the flux pinning. This fea- 
ture gives the new possibility for increasing the current- 
carrying capability of composite superconductors by op- 
timization the material structure without changing of its 
chemical composition. 
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